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Abstract 

Goal, Scope and Background. To strengthen the evaluative power 
of LCA, life cycle interpretation should be further developed. A 
previous contribution (Heijungs & Kleijn 2001) elaborated five 
examples of concrete methods within the subset of numerical ap¬ 
proaches towards interpretation. These methods were: contribu¬ 
tion analysis, perturbation analysis, uncertainty analysis, com¬ 
parative analysis, and discernibility analysis. Developments in 
software have enabled the possibility to apply the five example 
methods to explore the much-used ecoinvent'96 database. 
Discussion of Methods. The numerical approaches implemented 
in this study include contribution analysis, perturbation analysis, 
uncertainty analysis, comparative analysis, discernibility analysis 
and the newly developed key issue analysis. The data used comes 
from a very large process database: ecoinvent'96, containing 1163 
processes, 1181 economic flows and 571 environmental flows. 
Conclusions. Results are twofold: they serve as a benchmark to 
the usefulness and feasibility of these numerical approaches, and 
they shed light on the question of stability and structure in an 
often-used large system of interconnected processes. Most of 
the approaches perform quite well: computation time on a mod¬ 
erate PC is between a few seconds and a few minutes. Only 
Monte Carlo analyses may require much longer, but even then 
it appears that most questions can be answered within a few 
hours. Moreover, analytical expressions for error propagation 
are much faster than Monte Carlo analyses, while providing 
almost identical results. Despite the fact that many processes 
are connected to each other, leading to the possibility of a very 
unstable system and very sensitive coefficients, the overall re¬ 
sults show that most results are not extremely uncertain. There 
are, however, some exceptions to this positive message. 

Keywords: Contribution analysis; discernibility analysis; 
ecoinvent'96; key issue analysis; life cycle interpretation; per¬ 
turbation analysis; sensitivity analysis; uncertainty analysis 


1 Goal, Scope and Background 

Part of the ISO 14043 definition of life cycle interpretation 
is in terms of the analysis of results from inventory analysis 
and/or impact assessment (Anonymous 2000). To this aim, 
three elements are distinguished: 

• Identification of significant issues based on the results of the LCI 
and LCIA phases of LCA, 

• evaluation which considers completeness, sensitivity and consist¬ 
ency checks, 

• conclusions, recommendations and reporting. 

As LCA is a tool for evaluation and to support decisions, it is 
clear that interpretation is in some respect a crucial phase of 


LCA. It is therefore necessary that sound procedures and pow¬ 
erful methods are available to support the interpretation phase. 

Obviously, the approaches mentioned in ISO 14043 require a 
quite diverse type of activity, expertise, and method. Sensitivity 
checks belong to the domain of statistical processing, whereas 
reporting requires communicational skills. In a previous contri¬ 
bution (Heijungs & Kleijn 2001), we proposed to distinguish 
procedural and numerical approaches to life cycle interpreta¬ 
tion, and elaborated five examples of concrete methods within 
the subset of numerical approaches. These methods were: 

• contribution analysis, 

• perturbation analysis, 

• uncertainty analysis, 

• comparative analysis, 

• discernibility analysis. 

They were implemented in the educational software tool, 
Chain Management by Life Cycle Assessment (CMLCA 
< http://www.leidenuniv.nl/cml/ssp/software/cmlca >). and 
examples of their use for a fictitional case study were pre¬ 
sented. In the meantime, CMLCA has evolved in such a way 
that it is able to handle very large systems, and an import 
facility for the ecoinvent'96 database (3rd edition, Frisch- 
knecht et al. 1996) has been created. This means that the 
practicality of the five example methods can be submitted 
to a test for large systems and, at the same time, that the 
much-used ecoinvent'96 database can be submitted to a 
number of potentially useful approaches in life cycle inter¬ 
pretation. The present paper describes such an analysis. 

Meanwhile, the numerical approaches towards interpreta¬ 
tion have developed further as well. Heijungs & Suh (2002) 
introduce two more methods that can be classified as a way 
of interpretation: 

• key issue analysis, 

• structural analysis. 

These have also been implemented in CMLCA, and this pa¬ 
per also describes how the key issue analysis performs on the 
ecoinvent'96 system, and conversely how the ecoinvent'96 
system performs on this test. The structural analysis is left out 
of this paper because it is still in a too early stage of develop¬ 
ment. For similar reasons, another newly developed approach, 
structural path analysis (Treloar 1997, Suh 2003), which is a 
method to disentangle important pathways in a network of 
unit processes, has not yet been included in this paper. A de¬ 
velopment that has been included is that matrix perturbation 
theory has provided analytical methods for uncertainty analy¬ 
sis (Heijungs & Suh 2002, Heijungs 2002), which may re- 
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place the Monte Carlo-based numerical approaches. As we 
will see in a next section, an analytical method may well be 
superior to a numerical method with respect to speed of com¬ 
putation without being worse in other respects. 

This paper can be considered as a sequel to our earlier one 
(Heijungs & Kleijn 2001). We will largely follow the structure of 
that paper. As before, the CMLCA program can be downloaded 
to verify and expand on our results. Notice, however, that the 
data used is not provided by us, because they are the property of 
the publishers of ecoinvent'96. If you possess the original CD, 
entitled 'Okoinventare von Energiesystemen. 3. Edition', you can 
redo our analyses. We have translated the German names of proc¬ 
esses, flows and compartments into English. The original Ger¬ 
man names have been added for easy reference in the Appendix 
to this article. It should be mentioned that a new release of the 
ecoinvent database has become available, ecoinvent data vl.01, 
also referred to as ecoinvent 2000 < http://www.ecoinvent.ch >. 
The timing and its price determined that we used the 3rd edition 
of the ecoinvent'96 for our analysis. In principle, the procedures 
should work in the same way. Computation requirements may 
be higher due to the larger size of the matrices. 1 

Some of the methods for interpretation require a specification 
of the uncertainties of the data, such as the probability distribu¬ 
tion function with a characteristic, e.g. a normal distribution 
with a mean of 100 kg C0 2 and a standard deviation of 3 kg 
C0 2 . These data are not available for ecoinvent'96. As a solu¬ 
tion, we have, where required, introduced a normal distribu¬ 
tion 2 with the given number as mean and a standard deviation 
of 5% of that value (hence, a variation coefficient of 0.05). 

All methods can be used at different levels of analysis, viz. 
inventory analysis, characterisation, normalisation and 
weighting. To concentrate on the ecoinvent database, the 
discussion is restricted to the inventory analysis. The gener¬ 
alisation to higher levels of analysis is, however, straightfor¬ 
ward. Computation times may be longer, of course, and some 
of the results may be quite uncertain because the impact 
assessment data introduce additional uncertainties. 

An important part of the performance aspect is of course 
the computer time requirements. For each of the methods 


1 During the review and revision period of this paper, the ecoinvent vl .01 has 
become public, and an interface has been added to CMLCA. To give an impres¬ 
sion: solving a system of 2522 processes and flows requires a few minutes. 

2 Of course, the validity of the assumption of normality and the value of 5% may 
be discussed in itself. This is, however, outside the scope of the present paper. 


described, we will give approximate computation times, 
found on a computer that could now be described as old- 
fashioned: Pentium III, 667 MHz, 128 MB RAM. With an 
optimized state-of-the-art computer, it should be possible to 
find much smaller computation times. Before turning to the 
interpretation-oriented methods, we should say something 
on the preparatory work that is done by CMTCA. Depend¬ 
ing on the exact route, several intermediate results are cal¬ 
culated and stored for future reference. These are: 

• the inverse of the technology matrix, A- 1 (see Heijungs & Suh 
2002, p.17), of which the computation takes 120 seconds, 

• the intensity matrix A (see Heijungs & Suh 2002, p. 19), of which the 
computation takes about 140 seconds (provided A -1 is available), 

• two matrices with the variances of the process data. 

When these matrices have been calculated and are thereby 
available, most methods for interpretation run smoothly. All 
timing details given are based on the assumption that these 
intermediate results have already been calculated. 

Hereafter, we will discuss all methods for interpretation that we 
have tried. Each method is described in a separate section. The 
description of each method that has been discussed by Heijungs 
& Kleijn (2001) is divided into the following subsections: 

• introduction; review of the basic concept, 

• results for ecoinvent'96, 

• performance. 

The key issue analysis, which has not been discussed by 
Heijungs & Kleijn (2001), receives a more extensive descrip¬ 
tion. The paper concludes with a discussion of the results of 
the interpretation of the ecoinvent'96 database and a reflec¬ 
tion on the role and future of life cycle interpretation. 

2 Discussion of Methods 

2.1 Contribution analysis 

The contribution analysis decomposes the aggregated results of 
inventory analysis, characterisation, normalisation or weight¬ 
ing into a number of constituting elements. For the inventory 
analysis, this means that a certain inventory item, e.g. the sys¬ 
tem-wide carbon dioxide emission, is traced back to the share that 
the different unit processes in the system are responsible for. 

2.1.1 Results for ecoinvent'96 

With a reference flow of 1 TJ UCPTE electricity, a contribution 
analysis of atmospheric carbon dioxide is shown in Table 1. It 
is immediately clear that only a few processes are responsible 


Table 1: Contribution analysis for carbon dioxide (to air) ([E25] as shown in CMLCA) associated with the functional unit 1 TJ UCPTE electricity ([G181] 
Strom-Mix UCPTE). Contributions below 3% are not shown; this explains that the entries shown add up to 79% of the system-wide release only 


Process 

Value (kg) 

Contribution (%) 

[P522] lignite power plant (Germany) 

2.87E4 

21 

[P580] coal power plant (Germany) 

2.3E4 

17 

[P400] oil thermic electricity (Italy) 

1.28E4 

9 

[P509] electricity from gas power plant (Italy) 

9.06E3 

7 

[P581] coal power plant (Spain) 

8E3 

6 

[P526] lignite power plant (Greece) 

6.03E3 

4 

[P511] electricity from gas power plant (Netherlands) 

4.68E3 

3 

[P513] electricity from gas power plant (West Germany) 

4.43E3 

3 

[P583] coal power plant (France) 

4.47E3 

3 

[P584] coal power plant (Italy) 

3.69E3 

3 

[P585] coal power plant (Netherlands) 

3.65E3 

3 
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for a large part of this emission. In this case, the emissions from 
coal combustion in Germany are seen to contribute to almost 
40% of the system-wide CO, emission of UCPTE electricity. 

2.1.2 Performance 

With the aggregated inventory results available, a contribution 
analysis for ecoinvent is ready within less than one second. 

2.2 Perturbation analysis 

The perturbation analysis identifies sensitive parameters, i.e. 
input parameters of which a small change induces a large 
change in a selected result. For instance, one may try to find 
out a process for which data a small change in data will lead 
to a large change in the carbon dioxide emission. The factor 
that couples a small change in input to a change in output is 
referred to as the multiplier. Multipliers larger than 1 (or 
smaller than -1) indicate sensitive parameters, while multi¬ 
pliers close to 0 indicate insensitive parameters. The ration¬ 
ale for using a perturbation analysis instead or on top of an 
uncertainty analysis is that it allows the researcher to study 
inherent sensitivities, even for variables for which no uncer¬ 
tainty indication is known. 


2.2.1 Results for ecoinvent'96 

Again, a reference flow of 1 TJ UCPTE electricity was chosen, 
as well as the target flow of atmospheric carbon dioxide. Re¬ 
sults are shown in Table 2. Here, we see that most parameters 
of the system have a very small influence on the carbon diox¬ 
ide emission. There is only one sensitive parameter (directly 
related to the electricity mixing process), and almost all multi¬ 
pliers are between -0.25 and 0.25. There are two multipliers 
in the intermediate region, around ±0.4; these relate to the 
share of German electricity in the UCPTE mix. 

CMLCA also allows one to run the perturbation analysis 
for all environmental flows in one step. This yields a very 
long table, of which we have brought only a small part into 
this paper (Table 3). In Table 3, the first column indicates 
the environmental flow that is perturbed as the result of a 
deliberate perturbation of the coefficient relating to the proc¬ 
ess in column 2 and the flow (often economic, but some¬ 
times environmental) in column 3. One sees that waste heat 
(to water) has many large coefficients, some of which are 
even larger than 100. A multiplier of 30 means that a mod¬ 
est input uncertainty of 3% would propagate as an output 
uncertainty of about 100%, and a multiplier of 100 would 
propagate a 1% uncertainty as 100%. Suspended particles 
(to water) is the next sensitive substance, with multipliers 


Table 2: Perturbation analysis for carbon dioxide (to air) ([E25] in CMLCA) associated with the functional unit 1 TJ UCPTE electricity ([G181] in CMLCA). 


Multipliers between -0.1 and 0.1 (many) are not shown 


Process 

Economic/environmental flow 

Multiplier 

[P293] electricity (UCPTE mix) 

[G181] electricity (UCPTE mix) 

-1.02 

[P294] electricity (West Germany mix) 

[G182] electricity (West Germany) 

-0.434 

[P293] electricity (UCPTE mix) 

[G182] electricity (West Germany) 

0.433 

[P589] electricity from lignite power plant (Germany) 

[G236] electricity from lignite power plant (Germany) 

-0.213 

[P294] electricity (West Germany mix) 

[G236] electricity from lignite power plant (Germany) 

0.211 

[P522] lignite power plant (Germany) 

[G626] lignite power plant (Germany) 

-0.211 

[P589] electricity from lignite power plant (Germany) 

[G626] lignite power plant (Germany) 

0.211 

[P522] lignite power plant (Germany) 

[E25] carbon dioxide (to air) 

0.21 

[P285] electricity (Italy mix) 

[G173] electricity (Italy mix) 

-0.209 

[P293] electricity (UCPTE mix) 

[G173] electricity (Italy mix) 

0.209 

[P597] electricity from coal power plant (Germany 

[G414] electricity from coal power plant (Germany) 

-0.175 

[P294] electricity (West Germany mix) 

[G414] electricity from coal power plant (Germany) 

0.171 

[P580] coal power plant (Germany) 

[G666] coal power plant (Germany) 

-0.17 

[P597] electricity from coal power plant (Germany 

[G666] coal power plant (Germany) 

0.17 

[P580] coal power plant (Germany) 

[E25] carbon dioxide (to air) 

0.168 

[P400] oil thermic electricity (Italy) 

[G390] oil thermic electricity (Italy) 

-0.111 

[P285] electricity (Italy mix) 

[G390] oil thermic electricity (Italy) 

0.109 


Table 3: Perturbation analysis for all environmental flows associated with the functional unit 1 TJ UCPTE electricity ([G181] Strom-Mix UCPTE). Most 


multipliers between -100 and 100 are not shown 


Environmental flow 

Process 

Economic/environmental flow 

Multiplier 

[El42] waste heat (to water) 

[P690] water flow power plant (UCPTE) 

[G753] water flow power plant (UCPTE) 

-128 

[El42] waste heat (to water) 

[P690] water flow power plant (UCPTE) 

[El 42] waste heat (to water) 

128 

[El42] waste heat (to water) 

[P692] water barrage power plant (UCPTE) 

[G755] water barrage power plant 
(UCPTE) 

-110 

[El42] waste heat (to water) 

[P293] electricity (UCPTE mix) 

[G182] electricity (West Germany) 

-106 

[El42] waste heat (to water) 

[P294] electricity (West Germany mix) 

[G182] electricity (West Germany) 

106 

[El42] waste heat (to water) 

[P692] water barrage power plant (UCPTE) 

[El 42] waste heat (to water) 

105 


Many more multipliers for [El42] waste heat (to water) 


[E455] suspended particles (to water) 

[P218] decarbonized water 

[W340] residue decarbonization in storage 

-15.7 

[E455] suspended particles (to water) 

[P922] residue decarbonization in storage 

[W340] residue decarbonization in storage 

15.7 
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around 15. Thus, the really problematic sensitivities are easily 
identified with the perturbation analysis, even if the exact 
uncertainty details of the process data are unknown. Far 
much lower on the list (not shown in the Table 3) are ethyl¬ 
ene oxide to air, acenaphtene to water and acrilonitrile to 
water, with multipliers between +4 and +5. But, as we have 
seen in Table 2, there are also environmental flows (CO z in 
this case) that do not occur in this high sensitivity table and 
that have multipliers between -0.1 and 0.1. 

2.2.2 Performance 

Heijungs & Kleijn (2001) speculated that a numerical varia¬ 
tion as a way to address the perturbation analysis would be 
quite problematic for large systems. This is confirmed by our 
analysis. The ecoinvent data contains approximately 9000 non¬ 
zero entries, and a successive variation of each of these pa¬ 
rameters with a solution to the system requires some 300 hours. 
For the analytical approach (cf. Sakai & Yokoyama 2002), 
however, we are now much more optimistic. The perturba¬ 
tion analysis for C0 2 (see Table 2) requires not more than 20 
seconds, and an automated run for all environmental flows 
(see Table 3) is completed within 5 min. This means that within 
a few minutes one can have a full analysis of the sensitivity of 
the ecoinvent for a chosen reference flow. For a different refer¬ 
ence flow, the analysis must be carried out anew. 

2.3 Uncertainty analysis 

The uncertainty analysis is devoted to the systematic study 
of the propagation of input uncertainties into output uncer¬ 
tainties. If uncertainties of the process data are specified, for 
instance in the form of a Gaussian distribution with a cer¬ 
tain standard deviation that may differ per process data item, 
the uncertainty analysis will produce standard deviations or 
confidence intervals for the inventory results. 

There are two basic ways of running an uncertainty analysis: 
by random sampling and by analytical formulas for error propa¬ 
gation. A well-known form of random sampling is the Monte 
Carlo analysis, of which the basic procedure is as follows: 

• every input parameter is considered as a stochastic variable with 
a specified probability distribution, 

• the LCA-model is constructed with one particular realisation of 
every stochastic parameter, 

• the LCA-results are calculated with this particular realisation, 

• the previous two steps are repeated a number of times (the sam¬ 
ple size A/), 

• the sample of LCA-results is investigated as to its statistical proper¬ 
ties (such as the mean, the standard deviation, the confidence 
interval). 

There are various extensions to this basic setup of Monte Carlo 
analysis, known under names like Latin hypercube sampling 
and the Metropolis algorithm; see Morgan & Henrion (1990) 
and especially Liu (2003) for an extended discussion. 

The analytical approach starts from the idea that the influ¬ 
ence of perturbations may be approximated by differential 
calculus, and that the variances due to independent stochastic 
perturbations are additive. For instance, when a mathemati¬ 
cal relationship is specified as 

z = f(x,y) 


then, the variance of z is approximately 
Var( z> = (i^) var(x) + [j^j var(y) 

Explicit matrix equations for LCA enable a similar elabora¬ 
tion of the variance of LCA-results (see Heijungs & Suh 
2002, p. 140 ff). It is noteworthy (cf. Maurice et al. 2000) 
that the analytical approach does not require that the form 
of the probability distribution of the input parameters be 
specified. It suffices to specify the first two moments (mean 
and variance). Thus, the data requirements of an analytical 
approach are lower than those of a sampling approach. 

2.3.1 Results for ecoinvent'96 

Input uncertainties are not available for ecoinvent. We have 
assumed that all process data are characterised by a Gaussian 
distribution with a standard deviation that is 5% of the base¬ 
line value. Thus, if a certain input data item was specified as 
200 kg, we replace it by N(200, 10), meaning a normal dis¬ 
tribution with a mean of 200 kg and a standard deviation of 
10 kg. Of course, one may dispute the validity of this 
Gaussian assumption of 5%, especially as the default distri¬ 
bution for ecoinvent data vl.01 is lognormal (cf. Hofstetter 
1998). Such discussions are outside the scope of this paper, 
which concentrates on a mere demonstration of how as¬ 
sumed uncertainties propagate. 

With 100 Monte Carlo runs, the reference flow 1 TJ UCPTE 
electricity yielded results for carbon dioxide and some other 
flows as seen in Table 4. Here, baseline indicates the result 
without stochastic calculations, while mean refers to the 
average result of the 100 Monte Carlo results. The column 
labeled Variation gives the coefficient of variation, the 
dimensionless ratio between the sample standard deviation 
and the mean; here it has been expressed as a percentage. In 
general, baseline and mean value agree well, and the agree¬ 
ment will increase when the number of Monte Carlo runs 
increases. For those inventory items for which the coeffi¬ 
cient of variation is very large, there may be quite some dis¬ 
crepancy between baseline and mean. This is an artifact of a 
small sample size: one cannot expect to cover the parameter 
range well with 100 runs if it is very large. Observe that for 
some flows (like HFC-134a to air and styrene to air) the 
two results do not agree in sign. This is also a clear sign of 
numerical instabilities and/or exceptionally large variation. 

We see that, although all input parameters in a very large 
and interconnected system have an uncertainty of 5%, the 
C0 2 -results have an uncertainty that is still surprisingly small 
(7%). It is also clear that the quite modest output uncertain¬ 
ties for C0 2 are not exceptional. Indeed, almost all output 
uncertainties are to be found somewhere between 7 and 15%. 
There are, however, a few environmental flows which have 
a much larger uncertainty. This is the case for the param¬ 
eters that were identified as sensitive in the perturbation 
analysis, most notably waste heat to water. Output uncer¬ 
tainties can be as large as a few thousand %, which pro¬ 
vides serious doubts on the validity of decisions that are 
(partly) based on such highly sensitive flows. 
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Table 4: Uncertainty analysis (on the basis of 100 Monte Carlo runs and on the basis of analytical formulas for error propagation) for a selected number 
of environmental flows associated with the functional unit of 1 TJ UCPTE electricity ([G181] Strom-Mix UCPTE) 


Environmental flow 

Baseline 

Mean 

(Monte Carlo) 

Variation 
(Monte Carlo; %) 

Variation 
(analytical; %) 

Unit 

[E20] cadmium (to air) 

0.00293 

0.003 

8 

7 

kg 

[E21] methane (to air) 

253 

262 

12 

12 

kg 

[E22] cyanide (to air) 

0.00016 

0.000164 

11 

11 

kg 

[E23] cobalt (to air) 

0.0175 

0.018 

9 

8 

kg 

[E24] carbon monoxide (to air) 

33.4 

34.2 

7 

7 

kg 

[E25] carbon dioxide (to air) 

1.37E5 

1.39E5 

7 

7 

kg 

[E26] chromium (to air) 

0.0136 

0.0139 

7 

7 

kg 

[E27] copper (to air) 

0.0374 

0.0384 

7 

7 

kg 

[E28] ethane (to air) 

1.2 

1.22 

9 

9 

kg 

[E29] ethylene (to air) 

0.129 

0.133 

9 

9 

kg 

[E30] acetylene (to air) 

0.00217 

0.00224 

11 

10 

kg 

[E31] iron (to air) 

1.51 

1.53 

8 

8 

kg 

[E32] hydrogen sulfide (to air) 

0.133 

0.136 

12 

12 

kg 

[E33] mercury (to air) 

0.00807 

0.00823 

8 

7 

kg 

[El34] HFC-134a (to air) 

-1.17E-6 

8.87E-8 

5057 

15 

kg 

[El35] HCFC-22 (to air) 

0.000222 

0.000228 

14 

13 

kg 

[E136] styrene (to air) 

1.49E-10 

—5.55E-12 

2597 

22 

kg 

[El37] dioxin/TEQ (to air) 

5.43E3 

5.54E3 

8 

8 

ng 

[El38] carbon tetrachloride (to air) 

2.9E-5 

2.98E-5 

10 

9 

kg 

[El39] chloroform (to air) 

1.34E-6 

1.42E-6 

17 

13 

kg 

[El40] vinyl chloride (to air) 

8.28E-6 

8.64E-6 

15 

12 

kg 

[El41] xylene (to air) 

0.858 

0.878 

8 

8 

kg 

[El42] waste heat (to water) 

-0.000569 

-0.000684 

1771 

2052 

TJ 

[El43] acenaphthene (to water) 

-1.25E-16 

9.64E-17 

462 

41 

kg 

[El44] acenaphthylene (to water) 

0.00785 

0.008 

12 

12 

kg 

[El45] acrilonitrile (to water) 

-1.19E-13 

9.18E-14 

458 

41 

kg 

[El46] bis(2-ethylhexyl) phthalate (to water) 

3.02E-9 

3.14E-9 

13 

11 

kg 

[El47] BOD5 (to water) 

0.0766 

0.079 

9 

9 

kg 

[El48] butyl benzyl phthalate (to water) 

6.46E-14 

-1.64E-15 

3794 

22 

kg 

[E149] 1,1,1-trichlorethane (to water) 

9.22E-8 

9.7E-8 

13 

13 

kg 


Table 4 also shows results for the analytical approach. The 
standard deviation is the square root of the analytically de¬ 
termined variance. Observe that no separate mean value can 
be calculated here. The coefficient of variation is the ratio 
between standard deviation and baseline result. 

The agreement between sampling approach and analytical 
approach is quite good, although each series of Monte Carlo 
trials has a random character, and the analytical approach is a 
first-order approximation only. For inventory items with a 
large coefficient of variation the agreement is bad, although 
both approaches tell us that the variation is large. It is difficult 
to tell which of the two approaches produces more correct 
results in this case. The sampling approach may need an ex¬ 
cessively large number of Monte Carlo runs (or a better strat¬ 
egy to cover sparsely populated regions of sample space), while 
the analytical approach may need to take more that just the 
first two moments into account, as the first-order approxima¬ 
tion may no longer be sufficient. Using the baseline instead of 
the sample mean for determining the coefficient of variation 
in the sampling approach may already provide an improve¬ 
ment. There are a few flows (like the CFC 134a) for which the 
analytical method suggests a rather small uncertainty (15%), 
where the sampling approach yields a huge uncertainty (5057%). 


We speculate that one outlier of the Monte Carlo sampling 
might be responsible for this. 

2.3.2 Performance 

One Monte Carlo run for the ecoinvent requires approximately 
30 seconds. An experiment with 100 runs thus requires some¬ 
what less than one hour. This is feasible, but not very practi¬ 
cal. And the number of Monte Carlo runs is typically set higher, 
e.g. to 500 or 1000. With our old-fashioned computer and the 
probably suboptimal algorithm in CMLCA, this would re¬ 
quire a full working day. Increased computer performance and 
smarter algorithms may decrease the computer time substan¬ 
tially. Indeed, we redid the calculations on a 2.4 GHz PC with 
18 seconds per Monte Carlo run. 

Far more promising, however, is the analytical approach to 
uncertainty analysis. CMLCA is able to carry it out in a few 
minutes time. Faster computers and more efficient coding 
may reduce this to a minute or so. 

2.4 Comparative analysis 

The comparative analysis is nothing more than a systematic 
place to simultaneously list the LCA results for different 
product alternatives. 
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Table 5: Comparative analysis for carbon dioxide (to air) ([E25 in CMLCA) associated with several alternative national electricity scenarios, correspond¬ 
ing to functional unit 1 TJ electricity. The UCPTE alternative has been set to 1 


A 

B 

CH 

E 

Ex-Ju 

F 

GR 

1 

L 

NL 

P 

UCPTE 

W-D 

0.51 

0.88 

0.038 

1.01 

0.96 

0.24 

2.51 

1.56 

3.94 

1.44 

1.3 

1 

1.43 


2.4.1 Results for ecoinvent'96 

Table 5 shows a tabular representation of the comparative 
analysis of 1 TJ electricity according to several national char¬ 
acteristics (Austria, Belgium, Switzerland, Spain, former Yu¬ 
goslavia, France, Greece, Italy, Luxembourg, Netherlands, 
Portugal, UCPTE, and Western Germany) for C0 2 . We see 
that there are quite some differences between Switzerland and 
Luxembourg, but that the difference between the Netherlands 
and Western Germany is very small. Note, however, that this 
size of the difference does not say anything to the significance 
of these differences. The section on discernibility analysis be¬ 
low will address the issue of significance. 

2.4.2 Performance 

Once the inverse technology matrix is available, CMLCA 
calculates the comparative analysis within a few seconds. 

2.5 Discernibility analysis 

The discernibility analysis combines the comparative analy¬ 
sis and the uncertainty analysis. It is based on comparing 
the product alternatives for a large number of Monte Carlo 
runs. In contrast to the uncertainty analysis, no analytical 
approach is available, at least not now, and probably never. 
Heijungs & Kleijn (2001) proposed a pure ranking-based com¬ 
parative analysis (based on rank-order statistics), but it can in 
principle be applied to the proper values as well. Thus, we 
may compute the LCA-results for the different product alter¬ 
natives per Monte Carlo run, and apply a f-test comparison 
among each pair of product alternatives. Or we may deter¬ 
mine the rank order of each product alternative per Monte 
Carlo run, and count the frequency of pair-wise preference. 


2.5.1 Results for ecoinvent'96 

As in the uncertainty analysis, it was assumed that all proc¬ 
ess data is associated with a Gaussian uncertainty of with a 
standard deviation of 5%. A Monte Carlo analysis of 100 
runs gives results for C0 2 as in Table 6. We see, for instance, 
that Austrian electricity (A) is discernible from all other na¬ 
tional electricities, better than most, worse than Swiss (CH) 
and French (F) electricity on C0 2 . Spanish electricity (E), 
however, is far less discernible. Its C0 2 is not significantly 
different from that of Belgium (B), former Yugoslavia (Ex-Ju) 
and UCPTE. The reader should observe that the results for 
Belgium and Spain are not significantly different, whereas a 
purely non-stochastic comparison (see Table 5) yields quite 
some difference between the two (0.88 and 1.01, respectively). 
This is an illustration of the fact that a large difference is not 
always a significant difference. The difference between Portu¬ 
gal and Italy is on the same order of magnitude (1.3 versus 
1.56) while being significant at the 97%-level. 

Values and 95%-confidence intervals of the different alter¬ 
natives can be displayed in a simple graph (Fig. 1). Observe 
here that the third and sixth bar from the right, Portugal 
and Italy, show quite some overlap of confidence interval, 
while the run-by-run ranking shows that Portugal beats Italy 
97 out of 100 times. 

2.5.2 Performance 

As for the uncertainty analysis, one Monte Carlo run takes 
approximately 30 seconds, hence a good analysis requires 
several hours. However, in contrast to the uncertainty analy¬ 
sis, no time-saving analytical approach is available at this 
moment. 


Table 6: Discernibility analysis for carbon dioxide (to air) ([E25] in CMLCA) associated with several alternative national electricity scenarios, correspond¬ 
ing to functional unit 1 TJ electricity. All data uncertainties have been assumed to be drawn from a Gaussian distribution with a standard deviation of 5%. 
The number of Monte Carlo runs is 100. The number 0.83 at the intersection of row 5 (E) and column 3 (B), for instance, indicates that E had a higher C0 2 - 
emission than B in 83% of the Monte Carlo runs. The asterisks indicate the conventional significance statistic: * means that p<0.05, ** that p<0.01, and *** 
that pcO.001 


R>C 

A 

B 

CH 

E 

Ex-Ju 

F 

GR 

1 

L 

NL 

P 

UCPTE 

W-D 

A 

- 

o'" 

1 

o'” 

o'" 

1 

o'” 

o'" 

o'” 

o'" 

o'” 

o'” 

o'” 

B 

1 

- 

1 

0.17 

0.24 

1 

o'" 

o'" 

o'” 

o'" 

o'" 

0.12 

o'” 

CH 

o'" 

o'" 

- 

o'” 

o'” 

o'” 

o'" 

o'" 

o'” 

o'” 

o'” 

o'” 

o'” 

E 

1 

0.83 

1 

- 

0.64 

1 

o'" 

o'" 

o'” 

o'” 

0.03* 

0.45 

o'" 

Ex-Ju 

1 

0.76 

1 

0.36 

- 

1 

o'” 

o'" 

o'” 

o'" 

o'" 

0.34 

o'” 

F 

o'" 

o'" 

1 

o'” 

o'” 

- 

o'” 

o'” 

o'” 

o'" 

o'” 

o'” 

o'” 

GR 

1 

1 

1 

1 

1 

1 

- 

1 

o'” 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

o'" 

- 

o'” 

0.79 

0.97 

1 

0.8 

L 

1 

1 

1 

1 

1 

1 

1 

1 

- 

1 

1 

1 

1 

NL 

1 

1 

1 

1 

1 

1 

o'" 

0.21 

o'” 

- 

0.84 

1 

0.6 

P 

1 

1 

1 

0.97 

1 

1 

o'" 

0.03* 

o'” 

0.16 

- 

0.99 

0.23 

UCPTE 

1 

0.88 

1 

0.55 

0.66 

1 

o'" 

o'" 

o'” 

o'” 

0.01" 

- 

0*" 

W-D 

1 

1 

1 

1 

1 

1 

o'" 

0.2 

o'” 

0.4 

0.77 

1 

- 
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carbon dioxide (to air) in kg/TJ 



Fig. 1 : Visualization of the 95%-confidence intervals for carbon dioxide 
(to air) ([E25] in CMLCA) associated with several alternative national elec¬ 
tricity scenarios (from left to right Austria, Belgium, Switzerland, Spain, 
Former Yugoslavia, France, Greece, Italy, Luxembourg, Netherlands, 
Portugal, UCPTE, and Western Germany), corresponding to functional 
unit 1 TJ electricity. All data uncertainties have been assumed to be drawn 
from a Gaussian distribution with a standard deviation of 5%. The number 
of Monte Carlo runs is 100. 

2.6 Key issue analysis 

The key issue analysis is not included by Heijungs & Kleijn 
(2001). It is discussed in more detail than the other methods. 


Thus, one may express the relative contribution of the un¬ 
certainty in x to the uncertainty in z as wax(z)lwar(z) x , and 
the relative contribution of the uncertainty in y to the uncer¬ 
tainty in z as var(z)/var(z) y . 

The analytical formulas for LCA are provided by Heijungs 
& Suh (2002). Calculating the derivatives requires some 
straightforward calculus, and again leads to explicit formu¬ 
las, that may be implemented in software for LCA. 

2.6.2 Possibilities 

Like the uncertainty analysis, the key issue analysis can be 
performed at the levels of inventory analysis, characterisa¬ 
tion, normalisation and weighting. Of course, as the level of 
aggregation progresses, so does the number of input data 
items, and hence the number of uncertain input data items. 

Although the key issue analysis resembles the contribution 
analysis, the types of directions for decomposition differ. In 
a contribution analysis, one can ask what is the contribu¬ 
tion to the weighted index of S0 2 from electricity production 
by way of acidification. In a key issue analysis, one can ask 
what the contribution to the uncertainty of the weighted in¬ 
dex is of the uncertainty in the S0 2 emission of electricity pro¬ 
duction, or of all the S0 2 emissions, or of all the unit proc¬ 
esses, or of the SO, acidification factor, or of all the S0 2 
characterisation factors, or of all the acidification factors. 


2.6.1 Basic concept 

The key issue analysis has been introduced earlier (Heijungs 
1996) and reappears in Heijungs & Suh (2002). The basic 
idea is that not only LCA results can be decomposed in a 
contribution analysis, but that the uncertainty of LCA re¬ 
sults can be decomposed as well. This is especially useful in 
an iterative LCA setup, where results of the key issue analy¬ 
sis direct the researcher's focus on those data items for which 
more accurate data should be gathered with priority. Per¬ 
haps confusingly, this type of analysis is sometimes referred 
to as sensitivity analysis (Saltelli et al. 2001). 

The key issue analysis is based on the additive property of 
variances in the analytical approach towards the uncertainty 
analysis. If, according to the equation in Section 2.3, 


2.6.3 Tabular and graphical representation 

The tabular and graphical representations of the results of a 
key issue analysis is the same as that for a contribution analy¬ 
sis (cf. Heijungs & Kleijn 2001). 

2.6.4 Restrictions and warnings 

For the contribution analysis, there were problems in inter¬ 
preting negative contributions, for instance due to avoided 
processes and negative characterisation factors. For the key 
issue analysis, no such problems exist, because variances are 
always non-negative. 

2.6.5 Results for ecoinvent'96 


var(z) = 



var(x) + 


var(y) 


then one may decompose the resulting variance into two 
contributions: 

var(z) = var (z) x + var(z) y 

where 


var(z) x = 



var(.v) 


and 

varfzfj, 



var(y) 


Table 7 shows a key issue analysis for the uncertainty in the 
atmospheric carbon dioxide emission associated with 1 TJ 
UCPTE-electricity. Three input uncertainties are responsi¬ 
ble for 69% of the total uncertainty in CO,. Of these three, 
two relate to the same process. Hence, one may reduce the 
uncertainty of the total C0 2 emission a factor of 2 by a 
more careful specification of only one process. 

2.6.6 Performance 

The time needed for carrying out a key issue analysis is only 
slightly larger than that for carrying out an analytical uncer¬ 
tainty analysis. Once an uncertainty analysis has been com¬ 
puted, CMLCA produces results for a key issue analysis 
within a few seconds. When no uncertainty analysis was 
carried out prior to the key issue analysis, the computation 
of the ecoinvent case is completed within a few minutes. 
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Table 7: Key issue analysis for carbon dioxide (to air) ([E25 in CMLCA) associated with the functional unit of 1 TJ UCPTE electricity ([G181] Strom-Mix 
UCPTE). Contributions to uncertainty below 2% are not shown; this explains that the entries shown accumulate to 85% of the total uncertainty only 


Flow 

Processes 

Variance (kg 2 ) 

% 

[G181] electricity (UCPTE mix) 

[P293] electricity (UCPTE mix) 

4.83E7 

51 

[G182] electricity (West Germany) 

[P293] electricity (UCPTE mix) 

8.74E6 

9 

[G182] electricity (West Germany) 

[P294] electricity (West Germany mix) 

8.78E6 

9 

[G173] electricity (Italy mix) 

[P285] electricity (Italy mix) 

2.04E6 

2 

[G173] electricity (Italy mix) 

[P293] electricity (UCPTE mix) 

2.04E6 

2 

[G236] electricity from lignite power plant (Germany) 

[P294] electricity (West Germany mix) 

2.08E6 

2 

[G626] lignite power plant (Germany) 

[P522] lignite power plant (Germany) 

2.08E6 

2 

[E25] carbon dioxide (to air) 

[P522] lignite power plant (Germany) 

2.06E6 

2 

[G236] electricity from lignite power plant (Germany) 

[P589] electricity from lignite power plant (Germany) 

2.13E6 

2 

[G626] electricity from coal power plant (Germany) 

[P589] electricity from lignite power plant (Germany) 

2.08E6 

2 

[G414] electricity from coal power plant (Germany) 

[P597] electricity from coal power plant (Germany) 

1.43E6 

2 


3 Conclusions 

3.1 Conclusions on the ecoinvent'96 database 

One could sometimes feel the fears that the intrinsic sensi¬ 
tivity of such a complex system as ecoinvent is such that it is 
not possible to obtain meaningful results in LCA. These fears, 
however, were never substantiated, and neither were they 
refuted, so that they remained more an urban legend than 
an established fact. Without pretending to have provided a 
definite analysis for the ecoinvent, it is now possible to draw 
a few better-founded conclusions. 

A contribution analysis is straightforward to execute, and 
reveals the underlying dependencies in an understandable 
way. In the present case, C0 2 -emissions from UCPTE elec¬ 
tricity highly depend on German coal-powered electricity. 

Somewhat related to the topic of perturbation analysis is 
the concept of the condition number (Heijungs & Suh 2002). 
It is an overall measure of the instability of the system. Its 
logarithm indicates the amount of significant digits that may 
get lost (Thisted 1988). The condition number of the 
ecoinvent data matrix is dramatically large: about 10 14 . This 
means that 15 significant digits of input information may be 
lost in the process of matrix inversion. And as most data in 
the ecoinvent database cannot be expected to be stated in 
more than one or two digits, one could claim that results 
obtained with this database have no meaning whatsoever. 
We should, however, keep in mind that the condition number 
can be regarded as an extreme worst-case indicator (Heijungs 
& Suh 2002, p. 138, Heijungs 2002). A more refined pertur¬ 
bation analysis confirms this. For most environmental flows, 
all multipliers are well between -1 and 1. Nevertheless, there 
are some multipliers on the order of 100. This indicates that 
an uncertainty of 1% propagates as an uncertainty of 100%. 
In other words, a quite modest uncertainty of 1% makes the 
output of 100% uncertain. Note, however, that this applies to 
only very few environmental flows, at least in the case of 
UCPTE electricity that was considered, where, for instance, 
the C0 2 -emission is reasonably certain. 

When an uncertainty analysis is carried out with a data un¬ 
certainty of 5 % for all coefficients in the system, the output 
uncertainties for most environmental flows are between 7 
and 15%. Only for the flows with extremely high multipli¬ 


ers, we find uncertainties in output that well exceed 100%, 
and that may rise as high as several thousand percent. These 
extreme uncertainties are remarkable, but perhaps even more 
remarkable is the relatively low uncertainty of most envi¬ 
ronmental flows. The uncertainty of C0 2 , for instance, is 
only 7%. This picture emerges both for the analytical ap¬ 
proach (by error propagation formulas) and the sampling 
approach (by Monte Carlo analysis). 

An even more optimistic result is obtained with the 
discernibility analysis. Even when individual confidence inter¬ 
vals of several product alternatives show a large degree of 
overlap, a simultaneous analysis shows that the alternatives 
may still be significantly different and hence discernible. 

The key issue analysis provides an easy way to find those 
coefficients of which a reduction of uncertainty is most 
needed. Of course, its practicality must be investigated by 
experience in a setting that transcends the ambitions of the 
research reported here. 

3.2 Conclusions on life cycle interpretation 

A first conclusion is the overwhelming number of methods 
for life cycle interpretation. If the precise difference in mean¬ 
ing between a perturbation analysis and an uncertainty analy¬ 
sis is already difficult to comprehend, the addition of a key 
issue analysis and the distinction between an analytical and 
a sampling uncertainty analysis may probably stupefy some 
people that have had no specific training in mathematics or 
statistics. That is certainly a problem in communicating the 
results of a life cycle interpretation, but one should also ac¬ 
knowledge that many numerical approaches towards inter¬ 
pretation are for the exclusive use of the LCA-practitioner: 
it supports the process of reiteration with respect to data 
collection. In that respect, it is to be regretted that guidance 
documents that stress the importance of iterative LCA set¬ 
ups do not mention practically available tools like the per¬ 
turbation analysis. Several methods, most notably pertur¬ 
bation analysis and key issue analysis, will probably not be 
used for direct communication with the commissioner or 
the intended audience. Contribution analysis and compara¬ 
tive analysis may be used in that way, while uncertainty 
analysis and discernibility analysis may be presented in a 
slightly more friendly form. 
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Table 8: Overview of the suggested use of the different approaches towards life cycle interpretation under different conditions of availability of uncertainty 
information and for different purposes 


Purpose of analysis 

Uncertainty information available (or assumed) 

None 

Only standard deviation 

Standard deviation and distribution 

Focus in data collection 

1. Contribution analysis 

2. Perturbation analysis 

Key issue analysis 

- 

Decision-support 

Comparative analysis 

Uncertain analysis (analytical) 

1. Uncertain analysis (sampling) 

2. Discernibility analysis 


The practical application of numerical approaches to life 
cycle interpretation may be enhanced by the development 
of protocols for the successive use of these approaches in a 
context of iterative LCA-practice or decision-making. Pertur¬ 
bation analyses and key issue analyses may directly steer the 
data collection efforts, as it points out which data items should 
be collected with the highest precision and for which data 
items crude estimations may well be sufficient. On the other 
hand, comparative analyses combined with uncertainty analy¬ 
ses and discernibility analyses may provide information to a 
decision-maker, as these methods provide information on the 
ranking of alternative options and the robustness of this rank¬ 
ing under uncertainties in data. Table 8 provides guidance in 
the use of the different approaches in various situations. 

To many practitioners, sampling methods are more appealing 
than analytical approaches. But they have a clear disadvan¬ 
tage with respect to computation time. When one Monte Carlo 
run takes 30 seconds, 1000 runs require a working day. The 
analytical approach can reduce this to a few minutes, while 
the results are basically the same. Moreover, we see that a 
number of 100 Monte Carlo runs is sufficient for most pa¬ 
rameters. Clearly, the message of Sonneman et al. (2003), that 
'a full quantitative uncertainty assessment is, so far, too time 
consuming to be applied on each LCI as a default', needs some 
adjustment, given the availability of analytical methods and 
combined with the fact that these methods have a lower data 
demand (only the standard deviation suffices). 

More generally, we see that there is a choice between analyti¬ 
cal and sampling methods for an ordinary uncertainty analy¬ 
sis, while the key issue analysis is entirely based on an analyti¬ 
cal method and the discernibilty analysis on a sampling method. 

3.3 General conclusion 

Although this paper reports only a few of the analyses that 
can be carried out, it clearly shows that numerical approaches 
to life cycle assessment can well be applied to a large system 
of interconnected processes, and that the analysis of the 
ecoinvent database reveals some of the subtleties of this sys¬ 
tem. In our view, life cycle interpretation should become a 
standard practice of any LCA that pretends to be more than 
a screening LCA. Sophisticated use of approaches that re¬ 
side within the domain of interpretation should give us an 
understanding of the quality and robustness of the decision- 
support that is obtained with LCA. The knowledge on sen¬ 
sitivity and stability thus obtained should in its turn steer 
the collection of data with a sufficiently high quality. The 
perturbation analysis and the key issue analysis offer the 
possibility to discriminate between those data items for which 
a low or moderate data quality is sufficient, and those data 


items for which a high data quality is essential. Thus, sensi¬ 
tivity and robustness analyses facilitate an efficient data col¬ 
lection protocol. This once more stresses the iterative na¬ 
ture of the LCA process, where inventory analysis and 
interpretation are strongly interwoven. 
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Appendix 

This appendix lists our translations into English of the German names for processes, flows and compartments in ecoinvent'96 as far as they occur 
in this paper. Observe that no clear distinction between names of processes and their main outflow has been made in ecoinvent'96 


Processes 


Label 

German name 

Enqlish name 

[P2181 

Wasser entkarbonisiert 

decarbonized water 

IP2851 

Strom-Mix 1 

electricity (Italy mix) 

IP2931 

Strom-Mix UCPTE 

electricity (UCPTE mix) 

[P2941 

Strom-Mix W-D 

electricity (West Germany mix) 

[P4001 

Strom oelthermisch 1 

oil thermic electricity (Italy) 

[P509] 

Strom ab Brenngas-Kraftwerk 1 

electricity from gas power plant (Italy) 

[P511] 

Strom ab Brenngas-Kraftwerk NL 

electricity from gas power plant (Netherlands) 

[P5131 

Strom ab Brenngas-Kraftwerk W-D 

electricity from qas power plant (West Germany) 

[P5221 

Brk Kraftwerk in D 

liqnite power plant (Germany) 

IP5261 

Brk Kraftwerk in GR 

lignite power plant (Greece) 

fP5801 

Stk Kraftwerk in D 

coal power plant (Germany) 

IP5811 

Stk Kraftwerk in E 

coal power plant (Spain) 

IP5831 

Stk Kraftwerk in F 

coal power plant (France) 

[P5841 

Stk Kraftwerk in 1 

coal power plant (Italy) 

IP5851 

Stk Kraftwerk in NL 

coal power plant (Netherlands) 

[P5891 

Strom ab Brk-Kraftwerk in D 

electricity from lignite power plant (Germany) 

[P597] 

Strom ab Stk-Kraftwerk in D 

electricity from coal power plant (Germany) 

IP6901 

Laufwasserkraft UCPTE 

water flow power plant (UCPTE) 

[P6921 

Speicherkraft UCPTE 

water barraqe power plant (UCPTE) 

P922 

Rueckstaende Entkarbonisierung in Reststoffdeponie 

residue decarbonization in storage 


Economic flows 


Label 

German name 

English name 

[G1731 

Strom-Mix 1 

electricity (Italy mix) 

[G1811 

Strom-Mix UCPTE 

electricity (UCPTE mix) 

[G182] 

Strom-Mix W-D 

electricity (West Germany) 

[G236] 

Strom ab Brk-Kraftwerk in D 

electricity from lignite power plant (Germany) 

IG3901 

Strom oelthermisch 1 

oil thermic electricity (Italy) 

[W3401 

Rueckstaende Entkarbonisierung in Reststoffdeponie 

residue decarbonization in storaqe 

[G4141 

Strom ab Stk-Kraftwerk in D 

electricity from coal power plant (Germany) 

IG6261 

Brk Kraftwerk in D 

liqnite power plant (Germany) 

IG6661 

Stk Kraftwerk in D 

coal power plant (Germany) 

[G7531 

Laufwasserkraft UCPTE 

water flow power plant (UCPTE) 

[G755] 

Speicherkraft UCPTE 

water barrage power plant (UCPTE) 


Environmental flows 


Label 

German name 

Enqlish name 

[E20] 

Cd Cadmium[Emissionen Luftl 

cadmium (to air) 

[E21] 

CH 4 Methan[Emissionen Luftl 

methane (to air) 

[E22] 

CN CyanidefEmissionen Luftl 

cyanide (to air) 

[E23] 

Co Cobalt[Emissionen Luftl 

cobalt (to air) 

[E24] 

CO Kohlenmonoxid[Emissionen Luftl 

carbon monoxide (to air) 

[E25] 

C0 2 Kohlendioxid[Emissionen Luftl 

carbon dioxide (to air) 

TE261 

Cr Chrom[Emissionen Luftl 

chromium (to air) 

[E27] 

Cu Kupfer[Emissionen Luftl 

copper (to air) 

[E281 

Ethan[Emissionen Luftl 

ethane (to air) 

[E291 

Ethen[Emissionen Luftl 

ethylene (to air) 

[E30] 

Ethin[Emissionen Luftl 

acetylene (to air) 

[E31] 

Fe Eisen[Emissionen Luftl 

iron (to air) 

[E32] 

H?S Schwefelwasserstoff[Emissionen Luftl 

hydroqen sulfide (to air) 

IE331 

Hg Quecksilber[Emissionen Luftl 

mercury (to air) 

(E1341 

R134a FKW[Emissionen Luftl 

HFC-134a (to air) 

IE1351 

R22 FCKWfEmissionen Luftl 

HCFC-22 (to air) 

[El 36] 

Styrol[Emissionen Luftl 

styrene (to air) 

[El 371 

TCDD-Aequivalente[Emissionen Luftl 

dioxin/TEQ (to air) 

[El 381 

Tetrachlormethan[Emissionen Luftl 

carbon tetrachloride (to air) 

[El 391 

Trichlormethan (Chloroform)[Emissionen Luftl 

chloroform (to air) 

[El 401 

Vinyl Chlorid[Emissionen Luftl 

vinyl chloride (to air) 

[El 41] 

Xylole[Emissionen Luftl 

xylene (to air) 

[El 421 

Abwaerme in Wasser[Emissionen Wasserl 

waste heat (to water) 

[E1431 

Acenaphthene[Emissionen Wasserl 

acenaphthene (to water) 

[El 441 

Acenaphthylene[Emissionen Wasserl 

acenaphthylene (to water) 

[El 451 

Acrilonitrile[Emissionen Wasserl 

acrilonitrile (to water) 

[El 461 

bis(2-ethylhexyl) PhtalatfEmissionen Wasserl 

bis(2-ethylhexyl) phthalate (to water) 

[El 47] 

BSB5[Emissionen Wasserl 

BOD5 (to water) 

[El 481 

Butyl Benzyl PhtalatfEmissionen Wasserl 

butyl benzyl phthalate (to water) 

[El 491 

Chlor. 1,1,1-Trichlorethan[Emissionen Wasserl 

1,1,1 -trichlorethane (to water) 

[E455] 

Schwebestoffe[Emissionen Wasserl 

suspended particles (to water) 


Compartments 


Label 

German name 

English name 

[M4] 

Emissionen Luft 

emissions to air 

M5 

Emissionen Wasser 

emissions to water 
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